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Using the Ritus eingenfuriction method we compute corrections to the electron mass mo in the 
presence of a moderate magnetic field eB ~ mo- From this we obtain an expression for the anomalous 
magnetic moment near the critic field. For this we solved numerically the Schwinger-Dyson equations 
in the rainbow approximation including all Landau levels without make any assumption respect to 
the field strength. 

PACS numbers: 



INTRODUCTION 



particularly we follow the last reference. 



Radiative corrections to the electron self-energy due to 
its interaction with an external magnetic field have been 
widely considered in the literature [TJ|11[31I1[31IS1|71[H]. 
The exact form of how the ground state energy of an 
electron may be shifted from Eq in the presence of an 
magnetic field was a source of misleading conclusions [3J 
[HI Hn] • Assuming that this effects can be calculated by 
adding to the Dirac Hamiltonian a term A/icr • B, where 
A/i is the Schwingcr anomalous magnetic moment, it was 
founded that the ground-state energy of an electron is [5] 

null] 



Eq = niQC |1 



3hB 



At: itIqC^ 



(1) 



where toq is the electron mass without magnetic fields. 
For B ^ lO^^Gauss this formula implies that the ground- 
state energy of an electron is close to zero, a fact that 
could have dramatical astrophysical and cosmological 
consequences 0. However, it was soon realized that this 
result can only be true for magnetic field intensities that 
are weak compared with B^ = Tn^/e ~ 10^'^Gauss, be- 
cause Eq ([ij only represents the linear term in the power 
expansion in eB [S], which is not the dominant one for 
strong magnetic fields. In recent publications [8j [T^ dy- 



namical effects of a strong magnetic field B » m^/e 
were calculated. These approaches, show an enhance- 
ment of the dynamical contribution to the dynamical 
mass. They claim that this is a new effect that could 
have important consequences. We confirmed this result 
with an independent calculation and show that this effect 
is well described by previous calculations which are based 
on the constant mass approximation [6J. Furthermore, we 
carry out an analogous calculation for a weak magnetic 
field eB < ttiq, from this we obtain an expression for the 
anomalous magnetic moment near the critic field Be- Re- 
cently this issue has been revived from a non-perturbative 
point of view for strong fields |14J. In our approach 
we have included higher Landau levels to allow magnetic 
field intensities of the same order or less than the electron 
mass using similar techniques to those used in the con- 
text of the magnetic catalysis [T ^ [TO] [Tfl [T5 t [TO] ET] ■ 



II. THE MASS FUNCTION 

It has been shown [H] that the mass operator in the 
presence of an electromagnetic field can be written as a 
combination of the structures 

7^n^ , af^-F^, , (F^,n-)2 , 75i^p.F^'^ (2) 

which commute with the operator (7 • 11)^, where 

= td^ - eA^-, F^, = a^^r - d.^A°;\ F""" = 

^e^^"^'^ F\r, cTfj.i' = ^[jtJ.:lu] and A""' is the external vec- 
tor potential. We take A"^* = _B(0, — j//2, x/2, 0), such 
that the magnetic field is B = Bz. 

The equation relating the two-point fermion Green's 
function G{x, y) and the mass operator M{x, y) in coor- 
dinate space reads 

7 • T\{x)G{x, y) - [d^x'M{x, x')G{x', y) = S\x - y), (3) 



where the mass function in the rainbow approximation 
in coordinate space is [23] 

M{x,x') =mQ5^{x - a;')lixi 

-ze'rGix,x')YD(^J{x-x'), (4) 

where l4x4 is the 4x4 spinorial identity matrix. In 
this work, we are interested in exploring the behavior of 
the mass function for arbitrary values of the magnetic 
field strength, thus we work near the Landau gauge [23] 

~ 0) to compare with the well known results in the 
weak and strong magnetic field limit. 

Schwinger |23j was the first to obtain an exact ana- 
lytical expression for the fermion Green's function in the 
presence of a constant electromagnetic field of arbitrary 
strength. However, in the presence of a constant external 
field, the fermion asymptotic states are no longer free par- 
ticle states, but instead are described by eigenfunctions of 
the operator (7^H^)^. Thus, the alternative representa- 
tion of G(x, y) proposed by Ritus [22, is more convenient 
for our purposes, since there the mass operator is diago- 
nal. A matrix constructed out of these eigenfunctions is 



2 



used to rotate the SDE Q to momentum space, yielding 



tr[n(rtp)] (Al[p||,np] — mo) = — 2ie 

•-k ''p 2 ^ 



M[{p- q)i\,nk] 



(27r)4 q2 + jg (p-q)2 2eBnfe-X2[(p_^) 
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n 2 
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where qi_ = (0,gi,g2,0), q\\ = (90,0,0,^3), 7 = eB/2, 
= Ep —p^ — 2eBn and L™ are Laguerre functions and 



n(np) 



1, if 77,p 7^ 

A((T=1), if rip = 0, 



with A((t) = 1(1 + icr7^7'^). This factor comes from the 
normahzation [17j 



d^x^p{x)^p,{x) = 6n^n^,Ss^s^,S\p\\ -pj|)n(np), (6) 



and bring out a factor of two in the left-hand side of 
Eq. ([5| for Hp — 0. In vacuum, working near the Landau 
gauge, we know that the wave function renormalization 
equals one. We assume that this is the case also in the 
presence of a magnetic field of arbitrary strength. Ac- 
cordingly, to get Eq. (|5| we have worked with the ansatz 
S(p) ~ M{p)l4x4, where is defined by 



d'^xd'^x'^p{x)M{x, x')^p,{x') 
^Kn,,^s,s^,S'{p\\ -pj|)n(np)S(p). (7) 



This assumption is good when considering a small exter- 
nal momentum and Up — [TH (TBI [MJ |^. Since the 
energy only depends on the principal quantum number 
rifc, we expect that M((p—q) || , n^.) should be independent 
of Sfc. We also assume that M{{p — q)\\,nk) is a slowly 
varying function of Uk and thus make the approximation 



Miip- q)\\,nk) ^ Miip-q)\\,nk = 0). 



(8) 



Hereafter, we employ the more convenient notation 
M{ki\,nk = 0) = A^(fc||) for generic arguments of the 
mass function. With these considerations, the sum over 
Sk can be computed by using the identity 



4 = (_l)'^-'(i±i):L-' 



s+/' 



(9) 



and the Eq 5.11.5.1 in Ref. namely 



k=0 



l + A 
A 



xLl^_^{x + y + xA+^), 



(10) 



yielding 



tr[n(np)] {M{p) - mo) = -2re^ E / 



e "-1 



M{{p-q)\\) 



'^^ + {(9 - p)] - 2eB{k + [ak + l)/2)} - M^^ 
X |2 + (1 - 0(1 - '5.,.J^ - (1 - 0'^.,.. ^1 
x(-l)-™(-l)'=L,t™(|i)L-(^-")(|l), 



(11) 



where k — n^ — ^^^ti ^nd m = Up— '^''2^ ■ It is worth men- 
tioning that after summing over Sk, the resulting equa- 
tion is the same as Eq. (50) in Ref. [16^ when considering 
the case Uk = 0, which corresponds to the strong field 
limit. Under a Wick rotation we have 



{p - q)j - 2eB{k + ^) - A12 
-1 

{p - g)| + 2eB{k + ^) + 
2eB Jo 

With this result the sum over k can be performed also 



by resorting to Eq. (10) yielding, after carrying out the 
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sums over and Op, 



tr[n(?ip)] {M{p) - mo) = +2e' 



(2^)4 Q2 



+ ^ 2 - (1 - 0^} + 



(13) 



where Q = and the argument of the Laguerre func- 
tions is 4(5^ sinLn(-). For consistency with the assump- 
tion that the mass function is a slowly varying function of 
the principal quantum number Eq. ([8|, we take Up = 0. 
In this case, Eq. ( 13 1 gets simplifiea since Lg = 1 and 
L^i = 0, so we get 



M{p) =mo + 



d^Q M{l^-Q)\\ 



(2^)4 Q2 




(1-0- 



(14) 



This result corresponds with the Eq.(8) in the reference 
[?T] . Following [5T] we get for arbitrary gauge fixing pa- 



rameter the Eq. ( 15 I 



M{y) = mo + 



^•1 dz 



4 Jn (27r)- 



■M[z\ I dxx^ 



2(1 + x) - (1 - (^a: + (1 - xfz - y{l - xf + -yz{l - xf 

X e^(i--)r[0, y{\ - x)] + (1 - e)(l - x) (^^ - 1 + \z{\ - x) 

1 



^{y-z) + 



2{l + x)~{l-i)[x+-z\l-xf 



X e^(i-)r[0, z(l - x)] (1 - 0(1 - ^) (-^ + ^^(1 - ^)) 



(15) 



A1^[Q||]-»e 



2eB 



where z = Q"^, y = p^/i'y and A = z + 
1 and r(a;,?/) is the incomplete gamma function. The 
set of approximations leading to Eq. ( 15 1 are the same 



that in [2T which make no reference to the strength of 
the magnetic field, therefore this equation is valid for 
arbitrary magnetic field intensities. If we take, ^ = 1, 
in Eq. (15) it reproduces Eq. (11) in !2T] with mo — 0. 
For strong magnetic fields it will be useful to compare 
the solutions to this equation with the solutions tho the 



corresponding equation for the lowest Landau level (LLL) 



M{p\\ 




MAiq\\) 4 



e^r[o, 
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ql + MHq\\) ^'477 

(16) 

to calculate this expression we have used the same as- 
sumptions of the kernel softness respect to the momen- 
tum variables that in Eq. (15 1. We can find numeri- 



cal solutions to the integral equation Eq. (15 1 for several 



4 



magnetic field intensities making an appropriate choice 
of the cutoff in every range. For weak magnetic fields, 
eB < ttiq, we take as cutoff the electron mass given that 
it corresponds to the highest momentum scale and make 
^ = 0. For weak magnetic fields the dimensional re- 
duction is missing and thus it is important to take into 
account the ultraviolet divergences since we are in a simi- 
lar regime to the vacuum. This fact makes the numerical 
work harder and conceptually some care should be taken. 
There, contrary to the strong magnetic field case, the dif- 
ference between the mass function and the bare electron 
mass mo is meaningless because it does not tend to zero 
when the magnetic field is turn off, leading to to a ficti- 
tious dynamical mass. We can make an rough estimate 
of this contribution by consider the QED one loop con- 
tribution to the self-energy in the perturbative case. In 
the Landau gauge we get 



Sm — 



^-phys — mo 



In 

An 1 



-3 



mo, 
(17) 

which is at least two magnitude orders higher than the 
magnetic field contribution. To remove this additive con- 
stant is sufficient consider only differences between the 
mass function. Numerical results are shown in the fig. [T] 

-(m-m„s) 



mo 
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FIG. 1: Dynamical electron mass as a function of the mag- 
netic field in units of mg. The continuous line corresponds 
to the linear contribution Eq. |T|. Dots correspond to the 
numerical solution including ALL Eq. (151. Here rriref = 
m{eB 
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-"ml) 



2.4 X 10"'' mo 



It is noticeable that the relation, Eq. ([TJ, remains as the 
leading term even for field intensities close to B^. = in^/e. 
This is an important result because, despite that the 
old formula Eq. (jlj is very well established, to the best 
of our knowledge no validation of this result including 
higher Landau levels is known. Besides fig. [T] shows that 
the assumption eB « m? is no too restrictive. It is 
very well known that the anomalous magnetic moment 
of the electron A// = 4"^^^. can be obtained from the real 
part of the mass operator which depends of the magnetic 
field [H m [551 [571 UHl US] ■ In our approach we can 



read it from fig. [T] by use the formula A/i = \5m{B)/B\, 
where dm{B) is the magnetic correction to the electron 
mass. We can see in fig. |2] that the anomalous magnetic 
moment is independent with the magnetic field up to 
magnetic fields intensities near 0.1 Be in agreement with 
previous calculations which only include lowest Landau 
levels [5]. For magnetic fields near Be the anomalous 
magnetic field is highly suppressed and the same concept 
of magnetic moment in this range should be revisited. 
In the fig. |2] there is an small gap between the analyt- 
ical value and our numerical solution, this should only 
be a numerical error. For the strong regime we take 
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FIG. 2: Anomalous magnetic moment as a function of the 
magnetic field in units of rriQ. Dots correspond to the anoma- 
lous magnetic moment A/i taking into account all Landau 
levels. Horizontal dashed line corresponds to Schwinger clas- 
sical result 2iihi£Am = 



A — V eB and calculate the solutions to Eq. (|15| for 
the interval. 



< eB < 



which corresponds to 
magnetic fields in the phenomenological range . In this 



m— /Href 
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FIG. 3: Dynamical electron mass as a function of the mag- 
netic field in units of niQ. The continuous line corresponds 
to the analytical solution of Gusynin and Smilga. Gray and 
Black dots correspond to the numerical solution including 
ALL and the LLL respectively. Here m^ef = m(-B = mo/e) 
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regime we make the calculation including higher Landau 
levels and compare with the LLL solution. The solutions 
for eB/rriQ = 1 for ALL and LLL are mj-ef — 2.2 x IO^'^toq 
and TOrof = 1.4 X 10"'^ mo respectively, shown that high 
Landau levels represent an important contributions at 
this field intensities. However if, as in the weak mag- 
netic case, we only consider mass differences the results 
are indistinguishable as is shown in fig. [3j This lead us 
to conclude that for this regime of coupling, the main 
effect of higher Landau levels is to include vacuum con- 
tributions without changing significantly the dependence 
of the mass with the magnetic field. For higher magnetic 
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FIG. 4: Dynamical electron mass as a function of the mag- 
netic field in units of m§. The continuous line corresponds 
to the non-perturbative analytical solution of Gusynin and 
Smil ga. Black dots correspond to the numerical solution to 
Eq. ( 16 1 and the dashed line corresponds to perturbative so- 
lution in the Ref. [3]. Here mrd = m(B = mo/e) 

fields we recover the agreement between the lowest Lan- 
dau level solution and the complete solution. Given that 
there are not important contributions by including higher 
Landau levels for strong magnetic field, we calculate the 
enhancement of the electron mass and compare the so- 
lutions with the previous solutions founded in the litera- 
ture [6 . From the fig.|4]we can conclude that the dynami- 
cal effects are well described by the perturbative solutions 
which does not take into account the mass dependence 



with the momenta. We note the excellent agreement be- 
tween the analytical solution in reference [6] and our nu- 
merical solution. This result should not be surprising, 
since despite that in reference [5] they begin calculat- 
ing perturbative corrections, at the end they get a self- 
consistent equation for the dynamical mass and hence 
non-perturbative effects are expected in their approach. 
In summary: we have shown that the Eq. ( 15 1 can repro- 



duce perturbative as non-perturbative results for differ- 
ent ranges in the magnetic field intensity. This equation 
is easily solved with lower computational power, allow- 
ing to carry out calculations for magnetic fields relevant 
for the typical astrophysical conditions. To avoid mis- 
leading conclusions when dealing numerically with higher 
Landau levels, we have shown the importance of distin- 
guishing vacuum contributions from magnetic field con- 
tributions. Numerically this can be easily implemented 
by only consider differences between the mass function 
for different magnetic field intensities. Furthermore, in 
this work we have shown that the magnetic independence 
of the anomalous magnetic moment is consistent up to 
magnetic field intensities near lO^^Gauss. At same time, 
our technique allows to make a precise description of the 
way in which this assumption is broken in the magnetic 
field interval lO^'^ — lO^^Gauss. A natural extension of 
this work is to consider magnetic field intensities in the 
interval 10*^ — lO^Gauss. For this an improvement of 
our numerical techniques is needed. Due to highly pre- 
cise measurements of the anomalous magnetic moment 
non-perturbative effects are strongly bounded. If we ex- 
trapolate our calculations to magnetic field intensities of 
lO^Gauss, where some experiments currently take place, 
it is not completely clear that this contributions become 
less than one part in 10^^ which represents the nowadays 
precision in the anomalous magnetic moment. This de- 
serves a future investigation and can be an excellent test 
for non-perturbative techniques. 
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